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This  numerical  study  systematically  investigates  the  influence  of  initial  injection  conditions  of  reactants 
on  flame  characteristics  from  a  parallel  multi-jet  burner  in  a  laboratory-scale  furnace.  In  particular, 
varying  characteristics  from  visible  flame  to  invisible  Moderate  or  Intense  Low-oxygen  Dilution  (MILD) 
combustion  is  explored.  Different  parameters  examined  include  the  initial  separation  of  fuel  and  air 
streams  (S),  air  nozzle  diameter  (Da),  fuel  nozzle  diameter  (D/),  and  air  preheat  temperature  (Ta).  The 
present  simulations  agree  qualitatively  well  with  previous  measurements  reported  elsewhere  for  two 
reference  cases  investigated  by  experiment.  A  number  of  new  and  significant  findings  are  then  deduced 
from  the  simulations.  For  instance,  all  S,  Da  and  D/  are  found  to  play  significant  roles  in  achieving  a  proper 
confluence  location  of  air  and  fuel  jets  for  establishing  the  MILD  combustion.  Particularly,  varying  Da  is 
most  effective  for  controlling  the  combustion  characteristics.  It  is  also  found  that  the  stability  limits  of 
the  non-premixed  MILD  combustion  varies  with  different  combustor  systems  and  inlet  reactant  prop¬ 
erties.  Moreover,  for  the  first  time,  several  analytical  approximations  are  obtained  that  relate  the  flue-gas 
recirculation  rate  and  the  fuel-jet  penetration  to  Da,  Dp  S  and  also  reactant  properties. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

The  method  of  MILD  (moderate  or  intense  low-oxygen  dilution) 
combustion  has  been  recognized  as  one  of  the  most  successful 
combustion  technologies,  developed  in  past  two  decades,  which 
can  save  energy  and  simultaneously  reduce  emissions  of  NO*  (nitric 
oxide)  in  particular.  This  technology,  sometimes  also  termed  as 
‘HiTAC  or  HTAC  (high  temperature  air  combustion)  [1],  ‘FLOX’ 
(‘flameless  combustion’  or  ‘flameless  oxidation’)  [2],  has  been 
rapidly  transferred  from  laboratory  tests  to  industrial  applications, 
in  particular,  for  steel  reheating  and  steel  heat  treatment  furnaces. 
It  has  played  a  significant  role  in  the  mitigation  of  combustion¬ 
generated  pollutants  (particularly  NOx)  and  also  of  greenhouse 
gases  by  achieving  high  thermal  efficiency  in  metallic  industry. 
There  is  no  reason  that  the  technology  is  just  limited  within  the 
above  industry  field;  instead,  it  has  an  exceptionally  great  potential 
to  benefit  other  industries  in  future. 
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The  MILD  combustion  occurs  volumetrically  and  is  controlled  by 
strong  flue-gas  recirculation.  More  specifically,  fuel  is  slowly 
oxidized  in  an  environment  where  oxygen  is  highly  diluted  by 
recirculated  exhaust  gases  while  temperature  is  beyond  the  local 
auto-ignition  point  of  fuel.  Under  such  conditions,  traditional 
flames  cannot  sustain  and  are  always  blown-off  due  to  either  high 
jet  velocity  or  strong  internal  flue-gas  recirculation.  Chemical 
reactions  occur  in  a  distributed  zone  with  significantly  reduced 
peak  temperatures  [2,3].  As  a  consequence,  the  flame  is  hardly 
visible,  the  temperature  distribution  is  rather  uniform,  the  net 
radiation  flux  increases  by  as  much  as  30%,  and  NO*  emissions 
reduce  dramatically  [1-3]. 

The  MILD  combustion  has  been  a  favourite  subject  of  R&D 
(research  and  development)  in  the  last  two  decades,  see,  e.g.,  [4-7]. 
Weber  et  al.  [4,5]  performed  a  series  of  experiments  on  the  MILD 
combustion  of  gaseous,  liquid  and  solid  fuels.  Their  measurements 
were  conducted  inside  a  furnace  operating  with  highly  preheated 
air  regime  (1300  °C).  Their  data  for  NG  (natural  gas)  showed  that 
a  substantial  improvement  in  net  flux  of  the  thermal  radiation  can 
be  achieved  under  the  MILD  combustion.  Both  the  mixing  pattern 
and  intensity  have  significant  effects  on  the  overall  performance  of 
the  furnace  firing  with  NG,  specifically  on  the  thermal  efficiency 
part.  The  combustion  process  of  light  oil  is  very  similar  to  that  of 
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natural  gas,  with  invisible  flames.  However,  combustion  of  heavy 
fuel  oil  and  coal  is  significantly  different,  and  the  flames  are  always 
visible  in  their  tests.  Weber  et  al.  [5]  pointed  out  that  further 
research  is  still  needed  to  optimise  the  burner  designs  for  max¬ 
imising  the  recirculation  and  mixing  inside  the  furnace.  Note  that 
their  test  furnace  of  square  cross-section  is  characterized  by 
a  burner  with  a  central  (vitiated)  air  nozzle  and  two  small  off-axis 
fuel  injectors,  and  an  exhaust  outlet  opposite  to  the  burner  inlets. 
Experimental  investigations  of  the  MILD  combustion  using  this 
burner-furnace  configuration  have  been  performed  also  by  other 
investigators,  e.g.,  Rottier  et  al.  [8,9]. 

Szego  et  al.  [7]  tested  a  furnace  of  different  configuration  for 
MILD  combustion.  In  their  setup,  the  fuel  is  injected  through  four 
off-axis  holes  and  the  air  enters  through  a  central  nozzle,  while  four 
exhaust  ports  are  located  in  between.  These  investigators  reported 
measurements  of  temperature  and  flue-gas  composition  from 
conventional  and  MILD  combustions  using  NG  and  LPG  (liquefied 
petroleum  gas)  as  fuels.  They,  as  well  as  Kumar  et  al.  [  10],  found  that 
air  preheating  is  not  required  to  achieve  MILD  combustion.  Their 
investigation  [7]  on  the  stability  and  operational  condition  found 
that  diluting  the  fuel  with  an  inert  gas  helps  achieve  the  MILD 
condition.  This  is  because  the  dilution  shifts  the  stoichiometric 
contour  to  the  high  scalar  dissipation  region  which  suppresses 
flame  propagation  and  leads  to  a  distributed  reaction  further 
downstream  [11].  A  certain  threshold  of  the  fuel-air  momentum 
ratio  GflGa  (-0.006  for  their  test  system)  was  claimed  to  be 
necessary  for  the  MILD  combustion  to  occur.  This  momentum  rate 
was  considered  to  ensure  the  penetration  of  fuel  jets  to  a  region 
classified  as  the  oxidation  zone. 

Mi  et  al.  [12,13]  used  a  different  burner  configuration  at  the 
furnace  indicated  above  in  Refs.  [6,7]  to  investigate,  both  experi¬ 
mentally  and  numerically,  the  effects  of  the  air-fuel  injection 
momentum  rate  and  premixing  the  two  on  the  MILD  combustion.  A 
number  of  different  patterns  of  partially  and  fully  premixed  reac¬ 
tants  were  found  to  work  very  well  for  MILD  combustion.  Their 
numerical  study  suggested  that  there  is  a  critical  momentum  rate  of 
the  inlet  fuel-air  mixture  below  which  the  MILD  combustion  does 
not  occur.  In  the  MILD  regime,  both  the  inlet  fuel-air  mixedness  and 
momentum  rate  impose  insignificant  influence  on  exhaust  emis¬ 
sions  of  the  MILD  combustion.  It  is  worth  noting  that  their  fuel  and 
air  jets  issued  respectively  from  a  central  and  an  annulus  nozzle 
(non-premixed)  or  both  from  an  annulus  nozzle  (premixed).  Later, 
Li  et  al.  [14,15]  further  investigated  impacts  of  various  particular 
injection  conditions  on  the  characteristics  of  fully  premixed  MILD 
combustion  from  a  single  jet  burner  in  the  same  furnace.  The 
injection  conditions  include  the  area  of  the  nozzle,  equivalence 
ratio,  thermal  input,  and  initial  dilution  of  reactants.  They  found 
that  all  these  parameters  have  significant  influence  on  the  MILD 
combustion.  It  was  also  revealed  that  the  premixed  combustion  can 
occur  only  when  the  injection  Reynolds  number  exceeds  its  critical 
value.  Moreover,  very  low  emissions  of  NOx,  CO  and  H2  were 
measured  for  the  premixed  MILD  combustion  under  various 
conditions. 

Very  recently,  Schaffel-Mancini  et  al.  [16]  have  developed 
several  concepts  of  the  power  station  boiler  firing  pulverized-coal 
in  the  context  of  the  three  key  points  of  MILD  combustion:  exis¬ 
tence  of  an  intensive  in-furnace  recirculation,  homogeneity  of  both 
the  temperature  and  the  chemical  species  fields,  and  uniformity  of 
heat  fluxes.  In  order  to  determine  the  boiler  shape  and  its  dimen¬ 
sions,  these  authors  performed  CFD-based  numerical  simulations 
to  optimise  both  the  distance  between  burners  and  location  of 
the  burner  block.  It  was  demonstrated  that  the  momentum  of 
the  combustion  air  stream  is  an  essential  design  parameter  driving 
the  in-furnace  recirculation.  Also,  the  impact  of  the  combustion 
air  temperature  on  boiler  performance  was  found  less  critical 


providing  that  the  intensive  in-furnace  recirculation  has  been 
created.  Schaffel-Mancini  et  al.  [16]  finally  concluded  that  MILD 
technology  could  be  a  realizable,  efficient  and  clean  technology  for 
pulverized  coal  fired  boilers. 

In  the  context  of  the  previous  work  indicated  above,  the  present 
study  has  been  undertaken  to  systematically  investigate  the  influ¬ 
ence  of  injection  conditions  of  separate  methane  and  air  streams  on 
the  establishment  of  diffusion  MILD  combustion  in  a  furnace 
similar  to  that  of  Refs.  [4,5,8].  The  investigation  is  performed 
through  RANS  (Reynolds-averaged  Navier-Stokes  equations) 
modelling.  Several  injection  conditions  are  selected  that  include 
the  air  preheat  temperature  (Ta),  separation  distance  between  fuel 
and  air  exits  (S),  air  nozzle  diameter  (Da)  and  fuel  nozzle  diameter 
(D/).  The  main  objective  of  the  study  is  threefold: 

(1 )  To  qualitatively  verify  our  simulations  for  MILD  combustion  by 
RANS-modelling  with  previous  measurements  of  Rottier  et  al. 
[8,9]; 

(2)  To  explore  whether  the  geometric  parameters  S,  Da  and  D/  or 
the  initial  properties  of  reactants  are  important  for  the  occur¬ 
rence  of  MILD  combustion;  and 

(3)  To  identify  the  key  factor(s)  that  determine  the  MILD 
combustion. 


2.  Computation  details 

2.1.  Furnace  configuration  and  reference  cases  for  MILD  combustion 

The  present  study  simulates  various  combustion  cases  occurring 
in  a  laboratory-scale  furnace,  some  of  which  have  been  tested 
experimentally  in  a  previous  study  of  Rottier  et  al.  [8,9].  The  detail 
of  the  furnace  is  given  in  Refs.  [8,9]  and  here  only  a  brief  description 
is  provided  on  it.  Fig.  1  shows  a  schematic  of  the  furnace.  The 
combustion  chamber  is  square  in  cross-section  of  500  mm  by 


(Dia.=  Da) 


Fig.  1.  Schematic  of  the  furnace  and  burner  configuration,  simulated  presently  and 
measured  previously  in  Rottier  et  al.  [8].  (a)  furnace  and  burner  configuration;  (b) 
nozzle  arrangement.  Dimensions  are  in  mm. 
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500  mm  and  1000  mm  high,  made  with  refractory  materials  for 
a  maximum  wall  temperature  of  1400  K.  Several  openings  are 
set  along  in  staggered  rows  on  the  four  sides  of  the  combustion 
chamber.  These  openings  can  receive  an  optical  window  or  a  probe 
stand  for  detailed  measurements  in  the  furnace,  or  a  refractory 
block  to  preserve  thermal  confinement  and  measure  wall  temper¬ 
ature.  The  burner  consists  of  a  central  air  nozzle  (Da  =  25  mm)  and 
two  off-axis  fuel  pipes  (Df  =  3  mm),  with  the  initial  air-fuel  sepa¬ 
ration  S  =  50.7  mm.  The  exhaust  outlet  (centreed)  on  the  furnace 
top  is  square  in  cross-section  with  dimensions  of  50  mm  x  50  mm, 
according  to  Masson  [17].  Combustion  air  can  be  preheated  up  to 
Ta  =  838  K  by  an  electric  heater.  The  reference  case  for  Rottier  et  al. 
[8]  is  the  one  under  the  operation  conditions  of  firing  methane 
(CH4)  at  a  thermal  input  of  P  =  18.5  kW,  an  equivalence  ratio  of 
</>  =  0.85  and  an  air  preheat  temperature  of  Ta  =  838  K.  They  con¬ 
ducted  a  set  of  measurements  by  decreasing  the  preheat  temper¬ 
ature  from  838  K  to  293  K  when  keeping  constant  other 
parameters.  For  those  tests,  the  total  mass  flow  rates  of  air  and 
methane  can  be  estimated  to  be  rha  ~  7.5  x  1CT3  kg/s  and 
rhf  ~  3.7  x  10-4  kg/s.  Rottier  et  al.  [8]  reported  that,  for  their 
reference  case  (Ta  =  838  I<),  the  wall  temperature  was  homoge¬ 
neous  and  equal  to  1320  K,  and  the  flame  was  invisible,  so  in  the 
MILD  combustion  regime.  They  also  stated  that  the  MILD 
combustion  applies  for  all  their  test  cases. 

For  the  present  study,  the  above  cases  of  Ta  =  838  K  and  293  K 
are  set  as  reference  cases  to  verify  the  RANS  modelling. 

2.2.  Injection  conditions  of  air  and  fuel  for  simulations 

The  present  study  is  aimed  at  investigating  the  correlation  of 
MILD  combustion  with  the  distance  between  the  reactant  exits  (S), 
air  preheat  temperature  (Tfl),  air  nozzle  diameter  (Dfl),  and  fuel 
nozzle  diameter  (Df).  Table  1  lists  all  the  cases  of  present  investi¬ 
gation.  Note  that  numerical  tests  are  performed  for  varying  Ta,  Da,  Df 
or  S  in  the  first  column  at  the  listed  values  when  keeping  constant 
those  parameters  in  the  second  column  of  the  table. 

2.3.  Computational  conditions  and  models 

The  commercial  computational  software  package  Fluent  6.3  [18] 
is  used.  Full  three-dimensional  structured  grids  are  constructed  to 
have  small  orthogonality  deviations.  The  final  grids  for  half  model 
of  all  the  cases  listed  in  Table  1  are  determined  to  be  about 
400,000-600,000  cells  after  checking  the  independence  of  simu¬ 
lated  results  on  the  grid  number  for  the  reference  case;  note  that 
geometrical  symmetry  is  used  to  reduce  the  computational  time. 

Fuel,  air  and  flue  gases  are  all  assumed  to  obey  the  ideal  gas  law. 
Their  specific  heats  are  a  function  of  temperature  (piecewise- 
polynomial)  [18].  Inlet  and  outlet  boundary  conditions  of  the 
computational  domain  for  each  case  are  set  as  constant  velocities 
(corresponding  to  fuel  and  air  injections),  whose  magnitude  may 
vary  when  changing  Ta,  Da  and  Df  and  constant  static  pressure 
(1.0  atm).  Moreover,  a  uniform  temperature  distribution  is  set  as 
the  boundary  condition  of  the  entire  inner  wall  of  the  furnace.  The 


Table  1 

Injection  Conditions  of  Air  and  fuel  for  P  =  18.5  kW  and  <fi  =  0.85. 


Injection  conditions  varied 

Injection  conditions  fixed 

Ta  (K)  =  293a,  423,  573,  838a 

Da  (mm)  =  10,  15,  20,  25,  60,  80 
D/(mm)  =  2,  3,  5,  7,  9 

S  (mm)  =  20,  50,  100,  200 

Df=  3  mm,  Da  =  25  mm,  S  =  50.7  mm 
Df=  3  mm,  S  =  50.7  mm,  Ta  =  838  K 

Da  =  25  mm,  S  =  100  mm,  Ta  =  838  K 
Df=  5  mm,  Da  =  20  mm,  Ta  =  838  K 

a  These  cases  are  referred  to  as  reference  cases  for  the  present  study. 


uniform  temperature  taken  is  justified  because  the  wall  tempera¬ 
ture  measured  was  indeed  homogeneous  [8].  For  the  reference 
case,  the  wall  temperature  of  1320  K,  reported  in  Ref.  [8],  is 
employed  for  the  calculation. 

The  standard  k—e  model  with  the  standard  wall  function  is 
used  for  modelling  the  turbulent  flow.  The  EDC  (eddy  dissipation 
concept)  model  with  two-step  chemical  kinetic  mechanisms  for 
CH4  and  one-step  mechanisms  for  C2H6  are  applied.  A  preliminary 
study,  aimed  at  checking  the  suitability  of  the  global  chemical 
kinetic  mechanism  for  the  simulation,  is  carried  out  by  comparing 
the  differences  of  the  temperature  distributions  between  a  2-D 
model  with  the  global  chemical  kinetic  mechanisms  and  a  2-D 
model  with  a  detailed  GRI-3.0  mechanism.  (Similar  checks  were 
also  performed  in  our  previous  studies  [14,15].)  The  2-D  model  is 
represented  as  the  central  cross-section  of  the  furnace  (y  =  0,  xz- 
plane).  It  is  found  that  the  result  for  mean  temperature  of  the 
global  mechanism  model  is  similar  to  that  obtained  by  the 
detailed  GRI-3.0  mechanism,  e.g.,  with  the  maximum  difference 
between  the  predicted  temperatures  being  about  20  K.  For  this 
reason,  the  global  mechanisms  (two-step  chemical  kinetic 
mechanisms  for  CH4  and  one-step  mechanisms  for  C2H6)  are  used 
for  all  the  present  simulations  in  order  to  reduce  the  computa¬ 
tional  time  and  make  simulations  with  3-D  furnace  affordable.  To 
lessen  the  computational  cost  of  time  integration,  the  present 
modelling  adopts  the  in-situ  adaptive  tabulation  (ISAT)  model  of 
Pope  [19]. 

To  model  the  gas  emissivity,  the  DO  (discrete  ordinate)  radiation 
model  with  WSGGM  (weighted  sum  of  gray  gas  model)  is  applied 
for  radiation  [20].  For  DO  radiation  model,  theta  divisions  and  phi 
divisions  define  the  number  of  control  angles  used  to  discretize 
each  octant  of  the  angular  space,  and  the  both  parameters  are  set  to 
5.  Theta  Pixels  and  Phi  Pixels  are  used  to  control  the  pixelation  that 
accounts  for  any  control  volume  overhang,  and  both  parameters  are 
set  to  2.  The  SIMPLE  algorithm  is  used  for  pressure  velocity 
coupling.  Moreover,  the  NO*  is  computed  by  taking  equilibrium 
conditions  for  the  reaction  02^20,  see  Ref.  [21].  The  NOx  forma¬ 
tion  rate  is  given  by 

d[NO]/dt  =  2kNIf[0]eq[N2], 

where  /<Nif  =  1.8  x  1011  exp[-38370/T]  m3/kmols. 

A  second-order  discretization  scheme  is  used  to  solve  all  gov¬ 
erning  equations.  Convergence  is  obtained  when  the  residuals  are 
less  than  10~6  for  energy  and  DO  intensity  and  10-5  for  all  other 
variables.  The  outlet  temperature  and  velocity  are  monitored  and 
their  variations  within  1  I<  and  0.1  m/s,  respectively,  are  allowed  for 
convergence  of  the  solution.  Simulations  are  performed  with  a  16- 
core  CPU  and  26  GB  memory  workstation  (HP-Z800). 

3.  Definitions  of  fuel-jet  penetration  and  relative 
recirculation  rate 

3.2.  Fuel-jet  penetration  (Zj) 

For  the  non-premixed  or  diffusion  combustion,  fuel  and  oxidant 
streams  must  discharge  individually  into  furnace  from  their  own 
nozzles,  flow  downstream  separately  over  some  distance  prior  to 
their  confluence,  and  then  merge  as  a  single  jet,  see  Fig.  2.  The  fuel- 
jet  penetration,  denoted  by  Zj,  is  defined  as  the  vertical  distance 
from  the  fuel-nozzle  exit  to  the  location  at  which  the  (weak)  fuel  jet 
has  just  attached  to  the  (strong)  air  jet.  In  the  present  paper,  Zj  is 
estimated  to  be  the  minimum  value  of  z  at  which  only  one  peak 
occurs  at  the  centreline  in  the  mean  velocity  (vz)  profile  across  the 
xy  plane. 
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where  the  relative  rate 
of  recirculation  is: 


m 

Kv  = - <— 

ma  +  m  f 


m 


down 


ma  +  mf 


Fig.  2.  Definitions  of  the  fuel-jet  penetration  and  relative  rate  of  the  exhaust-gas 
recirculation. 


3.2.  Relative  recirculation  rate  (I<v) 


The  intense  dilution  of  reactants  prior  to  their  global  combus¬ 
tion  reactions  is  one  of  the  vital  requirements  for  achieving  overall 
MILD  combustion  in  the  downstream  region  atz  >  Zj.  The  degree  of 
low-oxygen  dilution  may  be  measured  by  the  relative  entrainment 
rate  of  the  exhaust  gas  into  the  reactant  stream(s),  first  defined  by 
Wiinning  &  Wiinning  [2]  as  the  recirculation  rate  (using  this  term 
throughout  the  paper  for  consistency),  viz., 


Kv  =  mr/ (ma  +  rhf'j ,  (1) 

where  mr  is  the  mass  flux  of  the  exhaust  gas  recirculated  and  then 
entrained  while  (rha  +  rhf)  denote  the  total  mass  flow  rate  of  the 
injecting  reactants  (i.e.,  air  and  fuel).  For  the  present  cases, 
according  to  the  mass  conservation  law,  at  each  cross  section  of  the 
chamber,  the  mass  flux  of  the  exhaust  gas  entrained,  rhr,  must  equal 
to  the  downward  mass  flux  of  the  recirculated  exhaust  gas,  rhdown, 
namely 


mr(z)  =  mdown(z) 


J  J  pvz(x,y)dxdy 


A(Z) 


with  A(z)  being  the  area  for  vz  <  0.  For  the  same  flux  of  reactants, 
a  higher  value  of  rhdown  reflects  a  higher  entrainment  rate  of  the 
exhaust  gas  and  thus  a  greater  recirculation  rate  or  a  higher  value  of 
Kv.  It  is  worth  noting  that  a  greater  value  of  I<v  or  mr  corresponds 
to  a  lower  value  of  the  local  averaged  oxygen  concentration, 
Cq2,  and  thus  a  higher  degree  of  low-oxygen  dilution,  since 
Co 2  ~rh02/(rh02  +  mN2  +  rhf  +  mr)  and  ma~m02  +  rhN2. 

Here  a  note  should  be  taken  on  the  dependence  on  Zj  of 
Kv{z  =  Zj)  or  simply  I<v(Zj).  As  defined,  separate  fuel  and  air  streams 
globally  meet  and  mix  each  other  at  z  >  Zj.  In  the  region  upstream 


from  this  location,  both  fuel  and  oxidant  are  gradually  diluted  by 
the  recirculation  and  entrainment  of  hot  combustion  products  into 
their  individual  streams.  Whether  or  not  the  MILD  combustion  can 
be  established  depends  on  the  degree  of  this  dilution.  Apparently, 
as  Zj  increases,  Kv(Zj )  becomes  greater  and  so  occurs  a  higher 
dilution  that  enhances  the  occurrence  of  MILD  combustion. 


4.  Results  and  discussion 

4.1.  Present  predictions  versus  previous  measurements  for  the 
reference  cases 

Rottier  et  al.  [8]  found  that  the  MILD  combustion  mode  occurs 
for  the  case  with  Ta  =  838  K  and  also  that  when  reducing  the 
combustion  air  to  a  temperature  of  Ta  =  293  K.  To  check  the  validity 
of  the  RANS  modelling,  the  present  predictions  are  compared  with 
their  measurements  in  Fig.  3(a)  and  (b).  These  two  figures  show  the 
predicted  distributions  of  the  mean  temperature  (T)  in  the  central 
xz  plane  (y  =  0)  and  those  of  the  measurements  [8]  for  Ta  =  838  K 
and  293  K,  respectively.  Rottier  et  al.  [8]  used  a  thermocouple  of  50- 
pm  bare  B-type  to  measure  the  mean  temperatures  point  by  point 
at  numerous  locations  through  several  windows.  So,  they  could  not 
provide  a  set  of  perfect  measurement  data  since  their  point- 
measurements  were  made  obviously  at  different  times  over 
a  certain  period  in  which  the  whole  operation  was  impossible  to 
keep  at  an  invariable  state;  besides,  the  limited  number  of 
measured  locations  should  be  another  negative  factor  for  plotting 
temperature  contours.  It  is  hence  understood  that  their  resulting 
temperature  maps  are  quite  scattering  and  definitely  with  signifi¬ 
cant  errors.  Also,  it  is  unfortunate  that  they  did  not  provide  any  XY- 
plots  of  their  measured  data,  so  that  we  cannot  make  better 
quantitative  comparisons  for  the  mean  temperature.  Nevertheless, 
taking  into  consideration  all  the  factors  that  influence  the 
measurement  accuracy,  a  careful  check  of  the  images  suggests  that 
the  calculated  and  measured  results  of  the  mean  temperature 
qualitatively  agree  reasonably  well  for  both  cases  with  Ta  =  838  K 
and  293  K. 

Specific  observations  then  should  be  made  here  from  Fig.  3(a) 
and  (b).  As  two  small  fuel  jets  proceed  downstream,  they  are 
tending  to  the  central  air  jet  and  simultaneously  mixing  the 
surrounding  flue  gas  where  the  oxygen  concentration  [O2]  ~  4%.  In 
the  low-oxygen  and  fuel-rich  mixing  layer,  oxidation  takes  place 
slowly,  with  a  small  amount  of  heat  releasing,  and  thus  tempera¬ 
ture  rises  slightly  (reflected  by  a  thin  yellow  or  faint-red  layer  in 
Fig.  3(a)  or  (b)).  The  slow  oxidation  consumes  almost  all  the  oxygen 
so  that  [02]  ~  0  in  the  mixing  layer  of  the  fuel  jet  (see  Fig.  6  below). 
Before  their  confluence  (z  <  Zj\  the  main  fuel  and  air  jets  all  have 
entrained  a  large  quantity  of  hot  combustion  products  and  hence 
been  well  diluted  (and  heated  up)  through  their  mixing  layers. 
When  they  globally  mix  each  other  at  z  >  Zj,  no  fast  oxidation  is 
possible  and,  instead,  slow  reactions  take  place  volumetrically, 
therefore  no  flame  front  being  formed.  Under  such  circumstances, 
the  overall  MILD  combustion  is  established,  regardless  of  Ta  =  293  K 
or  838  K.  This  was  confirmed  experimentally  by  Rottier  et  al.  [8], 
who  showed  that  the  flame  front  is  invisible  across  the  entire 
chamber  and  therefore  the  combustion  is  certainly  in  the  flameless 
MILD  mode  for  the  two  cases.  Nevertheless,  both  the  calculations 
and  the  measurements  (see  Fig.  3a, b)  consistently  reveal  that  pre¬ 
heating  the  combustion  air  from  Ta  =  293  I<  to  Ta  =  838  I<  not  only 
shifts  the  central  main  reaction  zone  towards  the  burner  exit  but 
also  reduces  the  volume  of  the  reaction  zone.  Besides,  the 
maximum  temperature  in  the  main  reaction  zone  is  about  200  K 
(400  K)  higher  than  the  mean  temperature  of  the  recirculated 
exhaust  gas  for  Ta  =  293  K  (838  I<). 
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Fig.  3.  Present  predictions  versus  previous  measurements  [8]  of  the  mean  temperature  in  the  central  xz  plane  for  the  reference  cases:  D/=  3  mm,  Da  =  25  mm,  S  =  50.7  mm, 
P  =  18.5  kW  and  0  =  0.85.  Air  preheat  temperature:  (a)  Ta  =  838  K;  (b)  Ta  =  293  K.  Note  that  the  two  images  of  the  measured  data  [8]  were  reproduced  directly  from  the  original 
paper  [8]  and  thus  that  the  two  contour  scales  used  in  [8]  for  Ta  =  838  K  and  293  K  are  also  selected  for  presenting  the  present  data.  Dimensions  are  in  mm. 


4.2.  Dependence  of  the  diffusion  combustion  on  injection 
parameters 

4.2 A.  Effect  of  varying  the  preheat  temperature  of  combustion  air 
(Ta) 

To  investigate  the  effect  of  Ta,  Fig.  4(a)— (d)  display  contours  of 
the  calculated  mean  temperature  (T)  in  the  central  xz  plane  (y  =  0), 
while  Fig.  5(a)— (d)  shows  those  of  the  calculated  O2  concentration, 


for  Ta  varying  from  293  K  to  838  K.  As  Ta  is  increased,  chemical 
reaction  between  fuel  and  ‘warmer’  air  should  become  faster  so 
that  the  central  main  reaction  zone  at  z  >  Zj  shifts  upstream 
(towards  the  burner)  and  also  decreases  in  size.  In  other  words,  an 
increase  in  Ta  should  cause  the  high  temperature  (>1500  K)  zone 
(yellow  to  red  colors)  moving  closer  to  the  burner  while  the 
maximum  temperature  should  increase  and  occur  upstream  —  as 
is  indeed  seen  in  Fig.  4  (and  also  Fig.  3).  This  is  also  reflected  by  the 


Temperature  (K): 
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Fig.  4.  Effect  of  air  preheat  temperature  Ta  on  mean  temperature  distributions  in  the  central  xz  plane  (upper)  and  the  xy  plane  where  Tmax  occurs  (lower)  for  the  reference  case: 
Df  =  3  mm,  Da  =  25  mm,  S  =  50.7  mm,  P  =  18.5  kW  and  0  =  0.85.  Dimensions  are  in  mm. 
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Fig.  5.  Contours  of  the  02  concentration  in  the  central  xz  plane  ( y  =  0)  at  different  preheat  temperatures  for  Df=  3  mm,  Da  =  25  mm,  S  =  50.7  mm,  P  =  18.5  kW  and  0  =  0.85. 
Dimensions  are  in  mm. 


region  of  near-zero  oxygen  concentration  (a  blue  zone  surrounded 
by  the  contour  of  [O2]  =  0.5%)  shown  in  Fig.  5.  These  observations 
agree  qualitatively  well  with  the  measured  results  of  Rottier  et  al. 
[8]  for  both  the  mean  temperature  (see  Fig.  3)  and  the  mean  OH 
concentration  (see  their  Fig.  5).  Their  Fig.  5  shows  that,  as  Ta  is 
increased,  the  concentration  of  OH  increases  while  the  high  OH 
region  moves  upstream;  correspondingly,  the  local  heat  release 
rate  grows  and  the  main  reaction  zone  shifts  towards  the  burner 
exit. 

422.  Effect  of  varying  the  air  nozzle  diameter  (Da) 

Figs.  6-9  respectively  display  contours  of  the  mean  temperature 
(T),  O2  (left)/CH4  (right)  concentrations,  CO2  (left)/CO  (right) 
concentrations  and  NOx  concentrations  in  the  central  xz  plane 
(y  =  0)  obtained  for  Da  varying  from  10  mm  to  80  mm.  Evidently,  as 


Da  is  decreased  from  25  mm  (reference  case)  to  10  mm,  the  CO2  and 
CO  concentration  gradients  in  the  reaction  zone  reduce  signifi¬ 
cantly,  particularly  the  mean  temperature  becomes  extremely 
uniform  and  no  high-temperature  region  at  all  appears  to  exist 
across  the  furnace,  thus  reflecting  better  performances  of  the  MILD 
combustion;  note  that  the  combustion  for  Da  =  25  mm  is  in  the 
MILD  regime  based  on  the  experimental  result  [8].  However,  the 
situation  is  totally  different  when  Da  is  increased  from  25  mm  to 
60  mm  and  80  mm:  the  fast  oxidation  process  and  thus  high  heat 
release  or  high  mean  temperature  occur  in  a  significantly  narrower 
zone.  For  these  cases,  the  maximum  mean  temperature  Tmax 
reaches  as  high  as  1974  K  and  2142  K  respectively,  and  so 
the  traditional  combustion  with  visible  flames  is  expected  to 
occur.  Interestingly,  Fig.  8  shows  that  the  high  CO  region  (>1.0%) 
grows  considerably  as  Da  is  increased  from  25  mm  to  80  mm. 


Temperature  |K]:  800  900  1000  1100  1200  1300  1400  1500  1600  1700  1800 


-250  0  250 

Da  [mm]:  10  15  20  25  60  80 

Fig.  6.  Effect  of  the  air  nozzle  diameter  on  mean  temperature  distributions  in  the  central  xz  plane  for  Df=  3  mm,  S  =  50.7  mm,  Ta  =  838  K,  P  =  18.5  kW  and  0  =  0.85.  Dimensions  are 
in  mm. 
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Fig.  7.  Effect  of  the  air  nozzle  diameter  on  02  and  CH4  concentration  distributions  in  the  central  xz  plane  for  D/=  3  mm,  S  =  50.7  mm,  Ta  =  838  K,  P  =  18.5  kW  and  0  =  0.85. 


By  comparison,  the  CO  concentration  is  nearly  uniform  at  [CO]  ~ 
0.3%-0.4%  across  the  furnace  for  Da  <  25  mm.  Also,  as  demon¬ 
strated  in  Fig.  9,  the  NO*  formation  is  nearly  zero  across  the  entire 
furnace  for  Da  <  25  mm  and  very  significant  in  the  high-T  regions 
for  Da  =  60  mm  and  80  mm;  particularly  interestingly,  the  exhaust 
emission  of  NOx  for  the  traditional  combustion  with  Da  =  80  mm  is 
about  300  times  that  for  the  MILD  combustion  even  with 
Da  =  25  mm  (not  Da  <  20  mm). 

Moreover,  a  careful  inspection  to  Fig.  6  finds  that  the  exhaust 
temperature  at  the  furnace  outlet,  Tou t,  takes  the  highest  value  for 
Da  —  80  mm  and  successively  lower  for  Da  =  60  mm,  25  mm  and 
20  mm.  For  Da  20  mm,  Tgyj  is  nearly  invariable.  This  variation 
results  perhaps  from  a  decreased  rate  of  heat  transfer,  out  of  the 
chamber,  caused  by  slower  recirculation  of  the  exhaust  gas  due  to 
increasing  Da.  In  addition,  the  present  calculations  find  higher 
emissions  of  the  exhaust  O2  and  CO  for  Da  >  25  mm  than  for 
Da  <  25  mm,  suggesting  that  the  conventional  combustion  may  be 
less  efficient  than  the  MILD  combustion. 

The  above  observations  may  be  explained  as  follows.  In  general, 
for  the  same  mass  flux  of  the  central  air  jet  (rha),  as  Da  is  increased, 
the  air  jet  injection  momentum  rate  Ga  decreases  (cf.  Eq.  (3)  in 
Section  5)  and  consequently  the  air  jet  entrainment  (then  the 
recirculation  of  exhaust  gases)  reduces  (which  will  be  quantified 
later).  This  obviously  does  not  benefit  for  the  establishment  of  MILD 
combustion.  There  is  another  negative  effect  of  increasing  Da : 
namely,  it  reduces  the  effective  separation  of  the  fuel  and  air  jets, 
i.e.,  Se  =  S  —  ( Da  +  Df)l 2,  and  thus  beneficial  for  the  traditional 
combustion  to  occur.  On  the  other  hand,  an  increase  in  Da  makes 
the  fuel  jets  to  become  relatively  ‘stronger’  and  thus  to  penetrate 
more  deeply,  see  Figs.  6-9.  However,  as  quantified  later  by  the 
recirculation  rate  (7<TV),  the  first  negative  effect  of  increasing  Da  is 
much  stronger  than  the  positive  one,  so  that  an  increase  in  Da 
actually  results  in  a  great  decrease  in  oxygen  dilution  by  the 


exhaust  gas  recirculated  prior  to  the  confluence  of  the  diluted  fuel 
and  oxygen  at  z  =  Zj,  thus  accelerating  the  reaction  rate  between 
CH4  and  O2. 

4.2.3.  Effect  of  varying  the  fuel  nozzle  diameter  (Df) 

Figs.  10  and  11  demonstrate  the  effect  of  varying  Df  on  the  mean 
temperature/oxygen  concentration  and  the  CH4/CO2  concentra¬ 
tions  in  the  central  xz  plane  for  the  cases  with  Da  =  25  mm  and 
S  =  100  mm.  Here  the  distance  between  the  air  and  fuel  jet  nozzle 
centres  is  selected  to  be  S  =  100  mm  and  Df  is  varied  from  2  mm  to 
9  mm  so  as  to  reduce  the  cross  influence  of  varying  Df  and  also  the 
effective  separation  of  the  fuel  and  air  jets,  i.e.,  Se  =  S  -  (Da  +  D/)/ 2; 
apparently  the  cross  influence  decreases  with  increasing  S  and 
decreasing  Da  or  Df.  As  Df  is  increased,  the  fuel  momentum  rate  G/ 
decreases  (see  Eq.  (3)  in  Section  5),  the  fuel  jets  weaken,  thus  their 
penetration  Zj  reduces  (indeed  evident  in  Fig.  10(a-e))  and, 
accordingly,  both  the  fuel  and  air  jets  are  less  diluted  by  the 
entrained  flue  gas.  In  other  words,  an  increase  in  Df  makes  both  the 
02  and  CH4  concentrations  to  become  locally  higher  so  that  CH4  is 
locally  oxidized  more  rapidly  in  the  region  upstream  of  Zj. 
Concurrently,  the  maximum  temperature  Tmax  due  to  combustion 
grows  and  thus  the  mean  temperature  distribution  becomes  less 
uniform.  Figs.  10  and  11  also  demonstrate  that  the  reaction  zone 
decreases  in  size  and  shifts  upstream  as  Df  is  increased.  Apparently 
from  Figs.  4  and  5  and  Figs.  10  and  11,  the  effect  of  increasing  Df  is 
roughly  similar  to  that  of  increasing  Ta;  namely,  both  drive  the 
combustion  to  depart  from  the  MILD  regime. 

4.2.4.  Effect  of  varying  the  separation  of  the  fuel  and  air  exits  (S) 

The  initial  separation  (S)  between  the  centres  of  the  fuel  and  air 

jets  should  be  another  important  factor  that  influences  the  MILD 
combustion  characteristics.  To  check  this,  Fig.  12(a)-(c)  present 
the  mean  temperature  and  oxygen  concentration  fields,  while 


Fig.  8.  Effect  of  the  air  nozzle  diameter  on  C02  and  CO  concentration  distributions  in  the  central  xz  plane  for  D/=  3  mm,  S  =  50.7  mm,  Ta  =  838  K,  P  =  18.5  kW  and  0  =  0.85. 
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Fig.  9.  Effect  of  the  air  nozzle  diameter  on  NOx  concentration  distributions  in  the  central  xz  plane  and  on  exhaust  NOx  emissions  for  Df  =  3  mm,  S  =  50.7  mm,  Ta  =  838  K, 
P  =  18.5  kW  and  0  =  0.85. 


Fig.  13(a)— (c)  show  the  CH4  and  CO2  concentrations,  over  the 
central  xz  plane  forS  [mm]  =  50, 100  and  200.  Obviously,  the  fuel  jet 
penetration  Zj  increases  with  increasing  the  separation  S.  It  follows 
that  increasing  S  causes  both  fuel  and  air  jets  to  entrain  more 
exhaust  gases  or  to  be  diluted  more  thoroughly  before  their 
confluence  at  z  =  Zj.  Hence,  as  S  is  increased,  not  only  the 
combustion  reaction  slows  down  but  also  the  reaction  zone  grows 
in  size,  so  that  the  temperature  becomes  more  uniform  and  its 
maximum  drops,  see  Figs.  12  and  13. 

4.3.  Effects  of  injection  parameters  on  Zj  and  I<v(Zj) 

Figs.  9-13  demonstrate  that  increasing  Df  and  S  can  make  the 
fuel  jet  to  penetrate  more  deeply,  before  its  confluence  with  the 
central  air  jet,  and  accordingly  the  main  reaction  to  occur  in 
a  lower-oxygen  diluted  zone.  It  is  deduced  that  a  sufficient  increase 
in  Df  or  S  and  thus  in  the  penetration  distance  Zj  may  enable  the 
combustion  regime  to  transfer  from  the  traditional  to  MILD  mode. 
However,  when  examining  the  effect  of  Da,  see  Figs.  6-8,  the 
increase  of  Zj  does  not  appear  to  be  always  beneficial  for  the 
establishment  of  MILD  combustion.  It  should  be  therefore  signifi¬ 
cant  to  quantify  the  influences  of  the  four  injection  parameters,  i.e., 
Da,  Df,  S  and  Ta,  on  both  Zj  and  I<v(Zj). 

Fig.  14(a— d)  shows  their  effects  on  Zj  for  P  =  18.5  kW  and 
0  =  0.85.  It  is  evident  that  Zj  basically  increases  as  Da  or  S  increases 
but  decreases  with  increasing  Df.  Also,  S  seems  to  have  the  strongest 
influence  on  Zj  while  Da  has  the  weakest.  Note,  nevertheless,  that 


increasing  Da  from  60  mm  to  80  mm  reduces  the  effective  sepa¬ 
ration  of  the  initial  fuel  and  air  jets,  i.e.,  Se  =  S  -  (Da  +  Df)/ 2,  from 
19.2  mm  to  9.2  mm,  thus  reducing  a  reduction  of  Zj,  see  Fig.  14(b), 
due  to  the  opposite  effects  of  Da  and  Se. 

Since  the  relative  rate  of  exhaust-gas  recirculation,  Kv(Zj), 
undoubtedly  depends  on  Zj,  it  thus  must  be  a  function  of  Da,  Df  and 
S.  This  is  indeed  confirmed  in  Fig.  15.  As  demonstrated,  I<v(Zj) 
decreases  significantly  with  increasing  Df  and  especially  Da  but 
increases  with  enlarging  S.  Apparently,  the  effect  of  varying  Da  on 
I<v(Zj)  and  thus  on  dilution  of  reactants  is  strongest,  noting  that 
I<v(Zj)  for  Da  =  10  mm  is  about  14  times  that  for  Da  =  80  mm.  In 
other  words,  varying  Da  from  10  mm  to  80  mm  (even  to  25  mm) 
change  the  combustion  structure  completely,  see  Fig.  6.  Moreover, 
the  effect  of  Ta  on  both  Zj  and  I<v(Zj)  is  found  insignificant 
(Fig.  14(d)).  This  is  consistent  with  the  finding  of  Schaffel-Mancini 
et  al.  [16]  that  the  impact  of  the  combustion  air  temperature  on 
the  performance  of  coal-firing  boiler  is  not  critical  if  the  intensive 
in-furnace  recirculation  has  been  created  for  MILD  combustion. 

Now  let  us  examine  the  effects  of  Zj  on  I<v(Zj)  and  Tmax,  the 
maximum  mean  temperature,  which  is  also  a  key  parameter  in 
MILD  combustion  as  it  is  responsible  for  NO*  production.  The 
results  obtained  from  variations  of  Da,  Df  and  S  are  presented  in 
Fig.  16.  To  make  a  meaningful  comparison  for  different  cases, 
normalization  of  Zj  by  Da  is  necessary  since  the  initial  momentum 
rate  of  the  central  air  jet  is  much  higher  than  that  of  the  fuel  jets 
and  thus  dominates  the  in-furnace  flow  structure.  It  is  evident  that, 
in  general,  as  Zj/Da  increases,  I<v(Zj)  increases  whereas  Tmax 


Fig.  10.  Effect  of  the  fuel  nozzle  diameter  Df  on  the  mean  temperature  and  02  concentration  in  the  central  xz  plane  for  Da  =  25  mm,  S  =  100  mm,  Ta  =  838  K,  P  =  18.5  kW  and 
0  =  0.85. 
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Fig.  11.  Effect  of  the  fuel  nozzle  diameter  Df  on  the  CH4  and  C02  concentrations  in  the  central  xz  plane  for  Da  =  25  mm,  S  =  100  mm,  Ta  =  838  K,  P  =  18.5  kW  and  0  =  0.85. 


decreases.  Namely,  the  low  value  of  rmax  corresponds  to  high  I<v(Zj) 
and  vice  versa.  However,  the  three  sets  of  data  for  either  Kv(Zj)  or 
Tm  ax  do  not  collapse  onto  a  single  curve.  In  particular,  the  data 
related  to  the  Da-variation  behave  quite  differently  when 
Da  <  20  mm.  As  Da  varies  from  10  mm  to  20  mm,  rmax  varies  little 
(1410  K-1490  I<)  while  I<v(Zj)  changes  significantly.  Close  looks  at 
Figs.  6-10  suggest  that  the  underlying  in-furnace  flow  and 
combustion  structures  for  these  three  cases  differ  substantially 
from  the  rest  cases  of  varying  Da. 


5.  Further  analysis  and  discussion 

Figs.  3-13  demonstrate  that  the  preheat  temperature  (Ta) 
and  the  geometric  parameters  (Dfl,  D/,  S)  of  air-fuel  injections  all 
have  significant  influence  on  the  diffusion  combustion  charac¬ 
teristics.  In  particular,  the  most  influential  factor  appears  to  be 
the  air  nozzle  diameter  (Dfl).  These  results  are  further  analyzed 
below. 

The  momentum  flux  (Ga)  of  the  central  air  jet  is  much  higher 
than  (~200  times)  that  of  each  fuel  jet  (G/)  so  that  the  air  jet  is 
expected  to  be  only  weakly  affected  by  the  fuel  jets.  These 
momentum  fluxes  can  be  expressed  as 


and  then  their  ratio  as 


Pa(lhl\2(Da\ 

Pf\maJ  \Df  J 


2 


7 


where  pa  and  pj  are  the  initial  air  and  fuel  densities,  ys  is  the  stoi¬ 
chiometric  air-fuel  ratio  (~  17.11  for  CH4),  and  Tf  is  the  initial  fuel 
temperature.  When  the  furnace  is  assumed  to  be  sufficiently  large, 
the  central  air  jet  may  be  ideally  regarded  as  a  single,  free  turbulent 
jet.  The  corresponding  fuel  jets  are  initially  parallel  to  and  eventually 
merge  to  the  air  jet  at  a  certain  downstream  point  (X/,  Zj).  That  is, 
each  fuel  jet  exhibits  a  curved  trajectory,  X(z),  towards  the  air  jet  and 
intercepts  the  mean  boundary  of  the  air  jet  at  z  =  Zj  (see  Fig.  2). 

We  consider  here  a  simplified  case  of  single  jet  where  the  effect 
of  buoyancy  is  negligible  and  the  jet  is  statistically  axisymmetrical. 
As  proved  well  by  previous  independent  investigations  (e.g. 
[22,23]),  the  mean  boundary  of  the  jet  can  be  described  by  the 
linear  relation 


B  =  Csz  (5) 

while  the  mass  flow  rate  is 

m  =  Ce(z-Zo)V0.25npa,Go;  (6) 


^PcJ^a 


and  Gf 


here  Cs  and  Ce  are  experimental  constants  (~0.2  and  0.32  for 
(3)  a  single  free  jet,  see  Ref.  [23]),  z0  is  the  z-location  of  the  virtual 
origin,  px  is  the  flue-gas  density,  and  G0  is  the  original  jet 
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momentum  flux.  Accordingly,  for  the  present  study,  the  ratio  of  the 
mass  entrainment  rate  (mr  =  m  -  m0 )  of  the  surrounding  gas  at 
z>Zj  to  the  original  injection  rate  of  reactants  (rh0  =  ma  +  rhf),  i.e., 
I<v  defined  in  Eq.  (1 ),  may  be  estimated  by 


Iilr  (Z-Z0)  1 

~ ~  ^e - - i 

ifi0  Da 


To  verify  Eq.  (7),  Fig.  17  presents  the  I<v  data  against  z/Da  over  the 
region  of  400  mm  <  z  <  800  mm  for  different  Da.  Surprisingly,  Eq. 
(7)  works  extremely  well  when  taking  Ce  ~  0.2,  for  the  combustor 
system  modelled,  despite  the  virtual  origin  z0  varying  with  Da.  It  is 
even  roughly  valid  for  Da  >  25  mm  (when  Ce  is  modified  slightly), 
where  the  furnace  should  not  be  considered  large  enough  for  the 
present  reactant  jets  to  be  ‘free’. 

In  this  context,  we  can  specify  the  analysis  of  Grandmaison  et  al. 
[24]  for  the  present  system  of  parallel  air  jet/fuel  jets  and  derive  the 
following  approximations  for  the  fuel  jet  penetration  Zj  and  the 
recirculation  rate  I<v(Zj): 


fG  \  1/4  s 

Kv{Zj)=C,CeU- )  Wa  +  C2,  (9) 

where  Cj  and  C2[=-(Cez0/Da+l)]  are  constants  that  should  be 
determined  by  experiment  or  numerical  simulation;  their  magni¬ 
tudes  should  depend  on  combustion  systems  of  investigation. 
Further,  substitutions  of  Eq.  (4)  into  Eqs.  (8)  and  (9)  attain  that 


Zj  *C10-1/2Ys“1/2  (Tf/Tay/4Dl/2D~V2S  (10) 

I<v(Zj)=C,CerV2v7y2(Tf/Tay/4D-y2Dfy2S  +  C2.  (11) 

Eqs.  (8)— (11)  can  be  checked  by  comparing  them  with  simulated 
results  for  different  geometric  parameters  Da,  Df,  and  S.  We  test  Eq. 
(8)  in  Fig.  18  using  the  data  obtained  under  varying  Da  and  Df.  It  is 
evident  that  this  approximation,  in  general,  works  well  with 
Ci  ~  4.9  no  matter  whether  the  momentum  ratio  Gf/Ga  is  varied  by 
changing  Da  or  Df.  There  is  an  exception  for  Da  =  80  mm  because  too 
large  a  dimension  of  the  air  jet,  confined  in  the  furnace,  has  strongly 
violated  the  conditions  for  Eqs.  (5)  and  (6).  Since  both  Eqs.  (7)  and 
(8)  hold  quite  well  (see  Figs.  17  and  18),  Eq.  (9)  should  be  valid 
approximately  as  well. 

In  this  context,  it  is  deduced  that  Eqs.  (10)  and  (11)  should  be 
applicable.  Indeed,  as  seen  in  Fig.  14,  Eq.  (10)  works  fairly  well  for  all 
the  geometric  parameters  while  Eq.  (11 )  is  good  for  S  but  not  so  for 
Df  and  Da,  especially  for  the  cases  with  Da  =  60  mm  and  80  mm.  In 
addition,  Eq.  (10)  provides  a  correct  trend  for  the  Ta  dependence  of 
Zj\  i.e.,  Zj  increases  with  decreasing  Ta,  which  is  consistent  with  the 
modelling  result. 

As  noted  early,  both  Kv  and  Zj  are  important  parameters  for  the 
occurrence  of  diffusion  MILD  combustion.  Increasing  them  yields 
a  higher  dilution  of  reactants  by  exhaust  gases  (e.g.,  C02,  H20,  N2) 
and  therefore  beneficial  for  establishing  the  MILD  combustion.  It  is 
thus  useful  for  burner  designers  to  find  key  factors  of  operation  that 
control  the  magnitudes  of  I<v  and  Zj.  Eqs.  (10)  and  (11)  suggest  that 
Zj  and  I<v  depend  not  only  on  burner  geometric  parameters  (e.g.,  Da, 
Df,  S)  but  also  on  the  initial  properties  of  reactants,  i.e.,  fuel  type, 
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Fig.  15.  Effects  of  Da,  Df,  and  S  on  Kv(Zj )  for  P  =  18.5  kW  and  0  =  0.85.  The  abscissa  is 
normalized  by  the  reference  quantity,  i.e.,  D*  =  25  mm,  Dj  =  3  mm,  or  S *  =  50.7  mm. 

air-fuel  ratio  (ys),  equivalence  ratio  (0),  and  initial  temperature  ratio 
(' Tf/Ta ).  However,  the  geometric  parameters  appear  to  be  more 
critical  for  controlling  Zj  and  I<v  and  thus  the  characteristics  of 
diffusion  MILD  combustion. 

Wiinning  &  Wiinning  [2]  provided  a  stability  diagram  for  non- 
premixed  MILD  combustion  of  methane  in  a  particular  furnace 
system  (Fig.  19).  By  relating  the  average  furnace  temperature  Tav  to 
I<v,  three  combustion  regimes  were  identified  and  termed  as 
Regions  A,  B  and  C  (as  indicated  on  the  plot).  Region  A  is  charac¬ 
terized  by  stable  conventional  flames,  Region  B  by  an  unstable 
transitional  regime  from  conventional  to  MILD  combustion,  and 
Region  C  by  the  MILD  regime.  Wiinning  &  Wiinning  [2]  estimated 
Kv  based  only  on  the  exhaust  gas  recirculated  into  reactants  before 
their  reaction  takes  place,  excluding  the  recirculation  of  hot  prod¬ 
ucts  into  the  reaction  zone  to  improve  combustion  stability.  This  is 
similar  to  the  present  calculation  of  I<v(Zj).  The  present  values  of  Tav, 
obtained  by  averaging  the  mean  temperatures  over  the  forward  or 


Fig.  17.  Flue-gas  recirculation  rate  (I<v)  at  different  z  for  Df  =  3  mm,  S  =  50.7  mm, 
Ta  =  838  K,  P  =  18.5  kW  and  0  =  0.85. 

positive  velocity  flow  region,  against  Kv(Zj )  are  also  shown  on  the 
same  diagram  of  Fig.  19.  Evidently,  only  the  case  for  Da  =  80  mm 
falls  in  Region  A  and  four  in  Region  B  while  the  rest  ten  are  located 
in  the  stable  MILD  region,  i.e.,  Region  C.  Our  modelling  results 
hence  suggest  that  the  detailed  classification  of  Wiinning  &  Wiin- 
ning  [2]  should  work  only  for  their  particular  case. 

The  above  suggestion  coincides  with  the  experimental  result  of 
Szego  et  al.  [7]  obtained  from  a  furnace  of  different  configuration. 
These  authors  found  that  the  lower  limit  for  stable  MILD  combus¬ 
tion  at  a  thermal  input  of  15  kW  is  G/Ga  ~  0.006,  irrespective  of 
firing  NG  or  LPG.  This  was  supported  by  their  thermal  field 
measurements.  Therefore,  they  believed  that  the  stability  limits  of 
Tav  -  I<v  depend  on  G//Ga.  Eq.  (9)  suggests  that  I<v(Zj)  °c  (G//Ga)1/4.  In 
other  words,  if  Eq.  (9)  is  valid  for  diffusion  combustion,  the  limit 
boundary  of  Tv  -  I<v{Zj)  for  the  region  “C”  should  depend  on  G//Ga. 
More  specifically,  it  varies  with  burner  geometric  parameters 
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Fig.  19.  Locally  averaged  mean  temperature  Tav  against  I<v  at  Zj  and  the  corresponding 
schematic  diagram  of  stability  limits  of  MILD  combustion  from  Ref.  [2]. 


(e.g.,  Da,  Df,  S),  and  also  fuel  type,  air-fuel  ratio  and  the  inlet 
temperature.  Szego  et  al.  [7]  established  stable  MILD  combustions 
in  their  furnace  with  Df  =  2  mm  (higher  Gf)  but  not  Df  =  3  mm 
(lower  Gf),  because  the  fuel-jet  penetration  for  the  latter  is  too  short 
to  enable  sufficient  dilution  of  reactants  for  MILD  combustion  to 
occur.  Hence,  the  stability  limit  for  MILD  combustion  varies  with 
different  non-premixed  combustor  systems  and  also  G/Ga. 


These  findings  should  also  apply  for  other  combustor  systems 
using  multi-jet  burners. 

(3)  Based  on  the  predicted  data,  the  analytical  approximations  Eqs. 
(8)  and  (9)  or  Eqs.  (10)  and  (11 )  appear  to  hold  reasonably  well 
for  varying  the  initial  fuel-air  momentum  ratio  (G/Ga)  or 
burner  geometric  parameters  (S,  Da  and  Df)  in  the  diffusion 
combustion  system  modelled.  Those  approximations  are 
believed  to  work  also  for  other  multi-jet  burner  systems  if 
taking  different  values  of  Ce,  Cj  and  C2.  If  they  are  indeed 
generally  valid,  even  qualitatively,  they  should  be  significant 
and  helpful  for  fundamental  research  and  particularly  tech¬ 
nology  development  of  the  MILD  combustion  from  a  multi-jet 
burner  system. 

(4)  The  critical  value  of  I<v(Zj )  for  the  stable  MILD  combustion  and 
thus  the  stability  limits  of  different  combustion  modes  vary 
with  particular  combustor,  burner  configurations,  the  injection 
momentum  ratio  G//Ga  and  reactant  properties.  Accordingly, 
the  schematic  diagram  of  Wiinning  &  Wiinning  [2]  obtained  for 
a  particular  combustion  system  is  not  suggested  to  be  used  as 
a  quantitative  benchmark  to  assess  whether  or  not  the  MILD 
combustion  truly  occurs  in  different  combustors  and/or  under 
different  reactant  injection  conditions. 
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6.  Conclusions 

A  systematic  study  by  the  RANS-modelling  has  been  carried  out 
on  the  diffusion  combustion  of  methane  and  air  issuing  from 
a  parallel  multi-jet  burner  at  a  laboratory-scale  furnace.  Different 
injection  conditions  have  been  considered  that  include  the  inlet 
separation  between  fuel  and  air  jets  (S),  air  nozzle  diameter  (Dfl), 
fuel  nozzle  diameter  (Df),  and  air  preheat  temperature  (Ta).  For 
those  cases  measured  by  Rottier  et  al.  [8],  the  modelling  predicts 
reasonably  well.  It  is  thus  believed  that  the  present  simulations 
should  work  qualitatively  for  other  cases  of  investigation.  More¬ 
over,  a  fundamental  analysis  performed  on  the  simulated  results 
has  led  to  several  interesting  relations. 

The  main  conclusions  drawn  from  the  present  study  are 
summarized  below: 

(1)  Two  parameters  are  found  important  for  investigating  the 
diffusion  MILD  combustion  of  a  multi-jet  burner:  i.e.,  the  fuel- 
jet  penetration  distance  Zj  and  the  relative  rate  of  recirculation 
of  exhaust  gases  occurring  at  the  fuel-oxidant  confluence  point, 
I<v(Zj).  Increasing  either  of  the  two  generally  yields  higher 
dilution  of  reactants  and  is  therefore  beneficial  for  establishing 
the  MILD  combustion. 

(2)  In  the  burner/furnace  system  modelled,  the  fuel-jet  penetra¬ 
tion  Zj  depend  significantly  on  Df,  Da  and  S.  Increasing  Df  and  S 
can  make  fuel  jets  to  delay  their  mixing  with  central  air  jet  and 
hence  the  main  reaction  to  occur  in  a  lower-oxygen  diluted 
zone.  However,  varying  Da  appears  to  be  the  most  effective  way 
to  control  the  magnitude  of  I<v  and  consequently  the  combus¬ 
tion  performance  and  stability  characteristics.  By  comparison, 
the  properties  of  reactants  have  weaker  effects  on  Zj  and  Kv(Zj). 
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